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Abstract 

We compare results of /^-factorization approach and next-to-leading order collinear-factorization 
approach for dijet correlations in proton-proton collisions at RHIC energies. We discuss correlations 
in azimuthal angle as well as correlations in two-dimensional space of transverse momenta of two 
jets. Some /cj-factorization subprocesses are included for the first time in the literature. Different 
unintegrated gluon/parton distributions are used in the /c^-factorization approach. The results 
depend on UGDF/UPDF used. For collinear NLO case the situation depends significantly on 
whether we consider correlations of any two jets or correlations of leading jets only. In the first case 
the 2^2 contributions associated with soft radiations summed up in the fej-factorization approach 
dominate at <ft ~ ir and at equal moduli of jet transverse momenta. The collinear NLO 2 — > 3 
contributions dominate over fcj-factorization cross section at small relative azimuthal angles as well 
as for asymmetric transverse momentum configurations. In the second case the NLO contributions 
vanish at small relative azimuthal angles and/or large jet transverse-momentum disbalance due to 
simple kinematical constraints. There are no such limitations for the /ct-factorization approach. 
All this makes the two approaches rather complementary. The role of several cuts is discussed and 
quantified. 
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I. INTRODUCTION 



The subject of jet correlations is interesting in the context of recent detailed studies 
of hadron-hadron correlations in nucleus-nucleus [lj and proton-proton j2j collisions. Those 
studies provide interesting information on the dynamics of nuclear and elementary collisions. 
Effects of geometrical jet structure were discussed recently in Ref.Q. No QCD calculation of 
parton radiation was performed up to now in this context. Before going into hadron-hadron 
correlations it seems indispensable to understand better correlations between jets due to the 
QCD radiation. In this paper we address the case of elementary hadronic collisions in order 
to avoid complicated and not yet well understood nuclear effects. Our analysis should be 
considered as a first step in order to understand the nuclear case in the future. We wish to 
address the problem how far one can simplify the calculation to be useful and handy in the 
nuclear case and yet realistic in the proton-proton case. 

In leading-order collinear-factorization approach jets are produced back-to-back. These 
leading-order jets are therefore not included into correlation function, although they con- 
tribute a big (~ ~) fraction to the inclusive cross section. The truly internal momentum 
distribution of partons in hadrons due to Fermi motion (usually neglected in the literature) 
and/or any soft emission would lead to a decorrelation from the simple kinematical con- 
figuration. In the fixed-order collinear approach only next-to-leading order terms lead to 
nonvanishing cross sections at <p ^ n and/or pi it ^ p2.t (moduli of transverse momenta of 
outgoing partons). In the /^-factorization approach, where transverse momenta of gluons 
entering the hard process are included explicitly, the decorrelations come naturally in a rel- 
atively easy to calculate way. In Figfj] we show diagrams illustrating the physics situation. 
The soft emissions, not explicit in our calculation, are hidden in model unintegrated gluon 
distribution functions (UGDF). In our calculation the last objects are assumed to be given 
and are taken from the literature. 




FIG. 1: Typical diagrams for fct-factorization approach to dijet production. 

The /c t -factorization was originally proposed for heavy quark production [4|. In recent 
years it was used to describe several high-energy processes, such as total cross section in 
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virtual photon - proton scatterin 
- heavy antiquark correlations [8 



production [14|, |15 



Higgs boson 



511, heavy quark inclusive production!!, B , heavy quark 
inclusive photon production (loL lllj ]. inclusive pion 
12| or gauge boson [13J production and dijet correlations 
in photoproduction [l6j] and hadroproduction (l7| . 

It is often claimed that the fc t -factorization approach includes implicitly some higher- 
order contributions of the standard collinear approach. This loose statement requires a 
better understanding and quantification. 

Here we wish to address the problem of the relation between both approaches. We shall 
identify the regions of the phase space where the hard 2^3 processes, not explicitly in- 
cluded in the leading-order /^-factorization approach, dominate over the 2^2 contributions 
calculated with UGDFs. We shall show how this depends on UGDFs used. 

We shall concentrate on the region of relatively semi-hard jets, i.e. on the region related 
to the recently measured hadron-hadron correlations at RHIC. Here the resummation effects 
may be expected to be important. The resummation physics is addressed in our case through 
the /ct-factorization approach. 



II. FORMALISM 

A. 2^2 contributions with unintegrated parton distributions 

It is known that at high energies, at midrapidities and not too large transverse momenta 
the jet production is dominated by (sub)processes initiated by gluons. In this paper we 
concentrate only on such processes. The region of forward/backward rapidities and/or pro- 
cesses with large rapidity gap between jets will be studied elsewhere. The cross section for 
the production of a pair of gluons or a pair of quark- antiquark can be written as 

da(hxh 2 -> jj) f d 2 k 1 jd 2 k 2 ,t 



d 2 Pi, t d 2 p 2 ,t 



Jdy ld y 2 d k ^ dk ^ ^L— \ M (gg ^ Jj)\ 2 
■ <5 2 ( k M + k 2 ,t - J> i,t - ~p 2,t)F(xi, klj)F(x 2 , kl t ) > ( 2 -l) 

where 

Xl = TMe+in + I^ e +, 2 (2.2) 

V s V s 

x 2 = —^e yi H j±e m . (2.3) 

The final partonic state is jj = gg, qq. If one makes the following replacement 

•FiOn.fcy - xi0i(xi)5(fcy (2.4) 

and 

•^2(^2, k 2 2 t ) -> x 2 g 2 (x 2 )5{kl t) (2.5) 

then one recovers the familiar standard collinear-factorization formula. 
The inclusive invariant cross section for g production can be written 

da{h x h 2 -»• j) f d 2 k u d 2 k 2± 

dyi d 2 Pu = 2 J dm —— ( - } k.=fe..+^ (2 - 6) 
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and equivalent ly as 



da(hxh 2 — ► j) f , d 2 ki >t d 2 k 2 j 



Let us return to the coincidence cross section. The integration with the Dirac delta 
function in (12.11) 

d yi dy 2 ^l^(...)S%..). (2.8) 

7T 7T 

can be performed by introducing the following new auxiliary variables: 

Qt = k l j4 + k 2 ,t , 

~q t = ~k u - ~t 2 + . (2.9) 
The jacobian of this transformation is: 

d(QuHt) _fi i\ (i n_ 2 . 2 = 4 . (210) 



l -l VI -l 



d{ k i, t , k 2 ,t) 
Then our initial cross section can be written as: 
do(h x h 2 -> QQ) 1 



,., j / d Vl dy 2 d 2 Q t d 2 q t (...) 6 2 (Q t - y M - ^ 2jt ) (2.11) 

d 2 pi,td 2 p 2 , t 4 J 

= i J d Vl dy 2 d\ (...) |^ =A = (2.12) 



= J / */icfr/2 'g*^£i (•••) |^ i= 7J t = (2.13) 

1 /" 1 * 

= - / dyidy 2 -dq 2 d<\) qt (...) |^ t= - . (2.14) 

Above P t = pi t t + p 2 ,t- Different representations of the cross section are possible. If one is 
interested in the distribution of the sum of transverse momenta of the outgoing quarks, then 
it is convenient to write 

d 2 pi,t d 2 p 2)t = ^d 2 P t d 2 p t = ^d<f) Pt P t dP t d<p pt p t dp t 

= ^2nP t dP t d(P pt p t dPt- (2.15) 
If one is interested in studying a two-dimensional map pxj x p 2> t then 

d 2 p u d 2 p 2:t = d(j) X Pi,tdp ljt d(f) 2 p 2 , t dp 2 , t • (2.16) 
Then the two-dimensional map in jets transverse momenta can be written as 



dv(Pi,uP 2 ,t) 
dpi >t dp2 it 



J dfadfc pi it p 2 , t J dy 1 dy 2 ^q t dq t d(f)q t (...) . (2.17) 



The integral over 0i and 02 must be the most external one. The integral above is formally 
a 6-dimensional one. It is convenient to make the following transformation of variables 



>1, VI, 



+ 2 , 0_ = 0i-0 2 ) , (2.18) 



where + G (0, 4tt) and 0_ G (— 2tt, 2tt). Now the new domain (0 + , 0_) is twice bigger than 
the original one (0i,02). The differential element 



>iacp 2 



The transformation jacobian is: 



> | u< > 2 \ 1 

~ 2 



(2.19) 



(2.20) 



Then 

d 2 p u d 2 p 2 , t = = Pi,tdpi, t P2,tdp 2 , 
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Pi,tdpi, t P2,tdp 2 ,t 2nd(p^ . (2.21) 



The integrals in Eq. (12. 17|) can be written equivalently as 



dpi tt dp 2 ,t 2 2 



y d<f>+d<f>- Pi,tP2,t J dy 1 dy 2 ^q t dq t d(f)g t (...) . (2.22) 



The first factor of \ comes from the jacobian of the transformation and the second \ is due 
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to the extra extension of the domain. 



as: 



By symmetry, there is no dependence on + and therefore the final result can be written 

d(T j Pl ' t ' P2 '^ = \ ■ \ ■ 4tt / rf0_ p u P2,t [ dy x dy 2 \q t dq t d<p qt (...) . (2.23) 
dpi,tdp 2 ,t 2 2 J J 4 

This 5-dimensional integral is now calculated for each point on the map pij x p 2>t . This 
formula can be also used to calculate a single particle spectrum of parton 1 and parton 2. 

The matrix elements for 2 — > 2 processes are discussed shortly in Appendix A. The 
analytical continuation of the standard on-shell matrix elements (see Appendix A) will be 
called in the following "on-shell approximation" for brevity. In Refs.[l3, 22| exact matrix 
elements for off-shell initial gluons were presented (see Appendix A). We have checked that 
the results obtained with the on-shell approximation and those obtained with the off-shell 
matrix elements are numerically almost identical. The deviations occur only for very virtual 
(large kt) gluons where the contribution to the cross section is small for majority of UGDFs. 

In the present calculation we shall include also components with gluon-quark and quark- 
gluon processes shown in FigfSJ In the next section we shall discuss how large are their 
contributions to the cross section. 



B. 2 — > 3 contributions in collinear- factorization approach 

Up to now we have considered only processes with two explicit hard partons. In this 
section we shall discuss processes with three explicit hard partons. In Figj3] we show a 
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FIG. 2: New ^-factorization contributions included in the present paper 




FIG. 3: A typical diagram for 2^3 contributions. The kinematical variables used are shown 
explicitly. 

typical 2 — > 3 process. We also show kinematical variables needed in the description of the 
process. We select the particle 1 and 2 as those which correlations are studied. This is only 
formal as all possible combinations are considered in real calculations. 

The cross section for hih 2 — > gggX can be calculated according to the standard parton 
model formula: 

Mhh, - m) = J**** ftfr.^fc.^ M99 - 999) (2.24) 
The elementary cross section can be written as 



dv(gg -> ggg) = — \M gg ^ ggg \ 2 dR 3 . (2.25) 
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The three-body phase space element is: 

j p d3 Pi d3 P2 d 3 p 3 4 4 

^ 3 = 2^(2^)3 2^(2^)3 2^(2^)3 (2?r) 5 ^ + ft "Pi "ft "ft) , (2-26) 
It can be written in an equivalent way in terms of parton rapidities 

dy x d 2 p u dy 2 d 2 p u dy 3 d 2 p u 44 
^ 3 = (4 7 r)(2 7 r)M4 7 r)(2 7 r)M4 7 r)(2 7 r) 2 (27r) 5 (Pa + P» " ^ " » " **) ■ (2-27) 

The last formula is useful for practical purposes. Now the cross section for hadronic collisions 
can be written in terms of 2 — > 3 matrix element as 

1 



da = dy x d p lit dy 2 d p 2 ,tdy 3 • ^CTT^p ^ ^1/1(^1, H f )x 2 f 2 (x 2 , fi f ) |-M 2 ^ 3 | 2 , (2.28) 
where the longitudinal momentum fractions are evaluated as 

Pl,i / . \ . P2,t , . \ | P3,t i . \ 

«i = ^exp(+yi) + — = exp(+y 2 ) + -^exp(+y 3 ) , 

V s V s V s 

x 2 = -^exp(-yi) + ^exp(-y 2 ) + exp(-y 3 ) . (2.29) 

V s V s V s 

Repeating similar steps as for 2 — > 2 processes we get finally: 



rfo- = 64 ^ 4 ^ 2 xxfiixx, /j, f )x 2 f 2 (x 2 , /j, f ) \M 2 ^ 3 \ 2 p u dp u p 2 , t dp 2 , t d(j)-dy 1 dy 2 dy 3 , (2.30) 

where 0_ is restricted to the interval (0, 7r). The last formula is very useful in calculating 
the cross section for particle 1 and particle 2 correlations. 



C. Unintegrated gluon distributions 

In general, there are no simple relations between unintegrated and integrated parton 
distributions. Some of UPDFs in the literature are obtained based on familiar collinear 
distributions, some are obtained by solving evolution equations, some are just modelled or 
some are even parametrized. A brief review of unintegrated gluon distributions (UGDFs) 
that will be used here can be found in Ref. |9[]. We shall not repeat all details concerning those 
UGDFs here. We shall discuss in more details only approaches which treat unintegrated 
quark/ ant iquark distributions. 

In some of approaches one imposes the following relation between the standard collinear 
distributions and UPDFs: 




(2.31) 



where a = xq or a = xg. 

Since familiar collinear distributions satisfy sum rules, one can define and test analogous 
sum rules for UPDFs. We shall discuss this issue in more detail in a separate section. 

Due to its simplicity the Gaussian smearing of initial transverse momenta is a good ref- 
erence for other approaches. It allows to study phenomenologically the role of transverse 
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momenta in several high-energy processes. We define a simple unintegrated parton distri- 
butions: 

F? auss (x, k 2 , /i 2 F ) = xpf l (x, n 2 F ) ■ f Gauss (k 2 ) , (2.32) 

where Pi° ll (x, ii 2 F ) are standard collinear (integrated) parton distribution (i = g,q,q) and 
fcaussik 2 ) is a Gaussian two-dimensional function: 



fGaussik 



J- I exp(-A ;i V2a 2 )I 



7T 



^2.33) 



The UPDFs defined by Eq.fl232D and fl2T3"3l is normalized such that: 

J F? auss (x, k 2 ,i? F ) dk 2 = xpf\x,iJ 2 F ) . 



(2.34) 



Kwiecihski has shown that the evolution equations for unintegrated parton distributions 
takes a particularly simple form in the variable conjugated to the parton transverse momen- 
tum. In the impact-parameter space the Kwiecihski equations takes the following relatively 
simple form 



dJ^Nsix^,^ 2 ) a s (fi 2 





27T/! 2 




- Fns(x 




_ a s (ix 2 ) 


dii 2 


2TT/J, 2 






dF G (x,b, n 2 ) 


_ a s (n 2 ) 


<9/! 2 


2tT/J 2 



dzPgg(z) 



Q(z - x) J ((l - z)pb) F NS ( -, b, ii' 



dz\ Q(z — x) Jq((1 — z)fib) 



P qq {z)F s l-,b,n 



x 



X 



- [ zP qq( Z ) + zP gq( Z )} ^s{x, b, fl 2 ] 

dz{Q{z-x)J Q {{l-z)ixb) P gq (z)F s (^,l 

- [zPgg(z) + ZPqg(z)) T G [X , 6, (/ 



(2.35) 



x 



+ P gg (z)J r G [-,b,ii 



We have introduced here the short-hand notation 

•F. ns = 3~ u — J~ui 3~d — J~& i 



(2.36) 



The unintegrated parton distributions in the impact factor representation are related to the 
familiar collinear distributions as follows 



F k (x, b = 0, ii 2 ) = -p k (x, n 2 ) . 



(2.37) 



On the other hand, the transverse momentum dependent UPDFs are related to the integrated 
parton distributions as 



xp k {x,/j, 



(2.38) 
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The two possible representations are interrelated via Fourier-Bessel transform 

poo 

T k (x,k 2 t ^) = \ dbbJ (k t b)F k (x,b,Li 2 ) , 



(2.39) 



F k {x, 6,/i 2 ) = / dk t k t J {k t b)T k (x,k 2 ^ 2 ) 



The index k above numerates either gluons (k=0), quarks (k> 0) or antiquarks (k< 0). 

While physically J- k (x, k 2 , /i 2 ) should be positive, there is no obvious reason for such a 
limitation for J- k {x, b, fi 2 ). 



In the following we use leading-order parton distributions from Ref.|34j as the initial 
condition for QCD evolution. The set of integro-differential equations in b-space was solved 
by the method based on the discretisation made with the help of the Chebyshev polynomials 
(see [23]). Then the unintegrated parton distributions were put on a grid in x, b and /i 2 and 
the grid was used in practical applications for Chebyshev interpolation. 

For the calculation of jet correlations here the parton distributions in momentum space 
are more useful. These calculation requires a time-consuming mult i- dimensional integration. 
An explicit calculation of the Kwiecihski UPDFs via Fourier transform for needed in the 
main calculation values of (xi, k 2 t ) and (x2, k 2 t ) (see next section) is not possible. Therefore 
auxiliary grids of the momentum-representation UPDFs are prepared before the actual cal- 
culation of the cross sections. These grids are then used via a two-dimensional interpolation 
in the spaces (x\, kf t ) and (x 2 , k 2 t ) associated with each of the two incoming partons. 



III. RESULTS 

Let us concentrate first on 2 — > 2 processes calculated with the inclusion of initial trans- 
verse momenta. We shall include the following four (sub)processes: 

• gluon+gluon — > gluon+gluon (called diagram A\, see Fig. la) 

• gluon+gluon — > quark+antiquark (called diagram A 2 , see Fig. lb) 

• gluon+(anti) quark — > gluon+(anti) quark (called diagram Bi, see Fig. 2a) 

• (anti)quark+gluon — > (anti)quark+gluon (called diagram B 2 , see Fig. 2b) 

Only first two were included recently in the k t - factorization approach [l7, 24]. The papers in 
the literature have been concentrated on large energies, i.e. on such cases when only gluons 
come into game. We shall show that at present subasymptotic energies (RHIC, Tevatron) 
also the last two must be included, even at midrapidities. Similar conclusion was drawn 
recently for inclusive pion distributions at RHIC [151 ] . 

In FigHJwe show two-dimensional maps in (pi,t,P2,t) for listed above subprocesses. Only 
very few approaches in the literature include both gluons and quarks and antiquarks. In the 
calculation above we have used Kwiecihski UPDFs with exponential nonperturbative form 
factor (b = 1 GeV -1 ) and the factorization scale ji 2 = {j> t ,min + Vt,max) 2 1 '4 = 100 GeV 2 . 

In FigJS] we show a fractional contributions (individual component to the sum of all four 
components) of the above four processes on the two-dimensional map (2/1,2/2). One point 
here requires a better clarification. Experimentally it is not possible to distinguish gluon 
and quark/ ant iquark jets. Therefore in our calculation of the (y\, 7/2) dependence one has to 
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FIG. 4: Two-dimensional distributions in p±j and p2 t t for different subprocesses gg — > (left 
upper) — > (right upper), (75 — > gg (left lower) and gg — > gg (right lower). In this calculation 
W = 200 GeV and Kwiecihski UPDFs with exponential nonperturbative form factor (60 = 1 
GeV -1 ) and n 2 = 100 GeV 2 were used. Here integration over full range of parton rapidities was 
made. 



10 




FIG. 5: Two-dimensional distributions of fractional contributions of different subprocesses as a 
function of y\ and yi for gg — > 55 (left upper) 55 — ► (right upper), gq — > 5(7 (left lower) and 
(751 — ► (751 (right lower). In this calculation W = 200 GeV and Kwiecihski UPDFs with exponential 
nonperturbative form factor and 60 = 1 GeV -1 were used. The integration is made for jets from 
the transverse momentum interval: 5 GeV < Pi,t,P2,t < 20 GeV. 
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symmetrize the cross section (not the amplitude) with respect to gluon - quark/ ant iquark 
exchange (yi — > 1/2, y 2 —* yi)- This can be done technically by exchanging t and u variables 
in the matrix element squared. While at midrapidities the contribution of diagram B\ + 
B 2 is comparable to the diagram A\, at larger rapidities the contributions of diagrams of 
the type B dominate. The contribution of diagram A 2 is relatively small in the whole phase 
space. When calculating the contributions of the diagram A\ and A 2 one has to be careful 
about collinear singularity which leads to a significant enhancement of the cross section at 
</>_=0 and y\ = y%, i.e. in the one jet case. This is particularly important for the matrix 
elements obtained by the naive analytic continuation from the formula for on-shell initial 
partons. The effect can be, however, easily eliminated with the jet-cone separation algorithm 
discussed in Appendix D. 




FIG. 6: The angular correlations for all four components: gg — > gg (solid), gg — ► qq (dashed) and 
9Q ~> 9Q = Q9 ~> Q9 (dash-dotted). The calculation is performed with the Kwieciriski UPDFs and 
60 = 1 GeV -1 . The integration is made for jets from the transverse momentum interval: 5 GeV 
< Pi,t,P2,t < 15 GeV and from the rapidity interval: -4 < 2/1,2/2 < 4. 

For completeness in FigJS] we show azimuthal angle dependence of the cross section for 
all four components. There is no sizeable difference in the shape of azimuthal distribution 
for different components. 

The Kwiecihski approach allows to separate the unknown perturbative effects incorpo- 
rated via nonperturbative form factors and the genuine effects of QCD evolution. The 
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Kwiecinski distributions have two external parameters: 

• the parameter 6 responsible for nonperturbative effects, such as primordial distribu- 
tion of partons in the nucleon, 

• the evolution scale fip responsible for the soft resummation effects. 

While the latter can be identified physically with characteristic kinematical quantities in the 
process fip ~ PitiPiti ^ ne ^ rs * one ^ s °f nonperturbative origin and cannot be calculated 
from first principles. The shapes of distributions depends, however, strongly on the value of 
the parameter 6 - This is demonstrated in Fig J7] for the gg — > gg subprocess. The smaller 
bo the bigger decorrelation in azimuthal angle can be observed. In FigJ7]we show also the 
role of the evolution scale in the Kwiecinski distributions. The QCD evolution embedded in 
the Kwiecinski evolution equations populate larger transverse momenta of partons entering 
the hard process. This significantly increases the initial (nonperturbative) decorrelation in 
azimuth. For transverse momenta of the order of ~ 10 GeV the effect of evolution is of the 
same order of magnitude as the effect due to nonperturbative physics. For larger scales of the 
order of \j? f ~ 100 GeV 2 , more adequate for jet production, the initial condition is of minor 
importance and the effect of decorrelation is dominated by the evolution. Asymptotically 
(infinite scales) there is no dependence on the initial condition provided reasonable initial 
conditions are taken. 

In FigJH] we show azimuthal-angle correlations for the dominant at midrapidity gg — > gg 
component for different UGDFs from the literature. Rather different results are obtained 
for different UGDFs. In principle, experimental results could select the "best" UGDF. We 
do not need to mention that such measurements are not easy at RHIC and rather hadron 
correlations are studied instead of jet correlations. 

Before we start presenting further more detailed results let us concentrate on NLO cal- 
culation 1 . In Figj9]we show the results of a naive calculation, on the {pi t t,P2,t) plane where 
soft divergences are shown explicitly. One clearly sees 3 sharp ridges: along x and y axes 
as well as along the diagonal. While the ridges along x and y axis can be easily eliminated 
by imposing cuts on p\ t and p2,t, i- e - on jets taken in the analysis of correlations. The 
elimination of the third ridge is more subtle and will be discussed somewhat later. Some- 
times asymmetric cuts on jet transverse momenta are imposed in order to avoid technical 
problems. 

Let us start from presenting the results on the plane (px,t,P2,t)- m Fig{10] we show the 
maps for different choices of UGDFs and for 2^3 processes in the broad range of transverse 
momenta 5 GeV < pi j t,P2,t < 20 GeV for the RHIC energy W = 200 GeV. In this calculation 
we have not imposed any particular cuts on rapidities. We have not imposed also any cut on 
the transverse momentum of the unobserved third jet in the case of 2 — > 3 calculation. The 
small transverse momenta of the third jet contribute to the sharp ridge along the diagonal 
Pi,t — P2,t- Naturally this is therefore very difficult to distinguish these three-parton states 
from standard two jet events. In principle, the ridge can be eliminated by imposing a cut on 
the transverse momentum of the third (unobserved) parton. There are also other methods 
to eliminate the ridge and underlying soft processes which will be discussed somewhat later. 

In FigJH]we show corresponding distributions in azimuthal angle 0_. Very different az- 
imuthal correlation functions are obtained for different UGDFs. The NLO azimuthal angle 
correlation function exceeds those obtained in the fc r factorization approach for 0_ < 90°. 



1 Please note that what we call here NLO, is called sometimes LO in the context of jet correlations 35]. 
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FIG. 7: The azimuthal correlations for the gg — * gg component obtained with the Kwiecihski 
UGDFs for different values of the nonperturbative parameter bo and for different evolution scales 
fi 2 = 10 (on line blue), 100 (on line red) GeV 2 . The initial distributions (without evolution) are 
shown for reference by black lines. 



When calculating dijet correlations in the standard NLO (2 — ► 3) approach we have taken 
all possible dijet combinations. This is different from what is usually taken in experiments 
3^ 1 , where correlation between leading jets are studied. In our notation this means p 3jt < pi t t 
and ps )t < p2,t- When imposing such extra condition on our NLO calculation we get the 
dash-dotted curve in FigJHJ In this case da/dtp^ = for 0_ < |7r. This vanishing of the 
cross section is of purely kinematical origin. Since in the /c t -factorization calculation only 
two jets are explicit, there is no such an effect in this case. This means that the region of 
0_ < |7r should be useful to test models of UGDFs. For completeness in FigJT2l we show 
a two-dimensional plot (pi t t,P2,t) with imposing the leading-jet condition. Surprisingly the 
leading-jet condition removes a big part of the two-dimensional space. In particular, regions 
with p 2 ,t > 2pi,t (NLO-forbidden regionl) and p\^ > %P2.t (NLO-forbidden region2) cannot 
be populated via 2^3 subprocess 2 . There are no such limitations for 2 — > 4, 2 — >• 5 and 
even higher-order processes. Therefore measurements in "NLO-forbidden" regions of the 



In LO collinear approach the whole plane, except of the diagonal pn = p2.t, is forbidden. 
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FIG. 8: The azimuthal correlations for the gg — > gg component obtained for different UGDFs 
from the literature. The Kwiecinski distribution is for 6q = 1 GeV -1 and fi 2 = 100 GeV 2 . 



(pi,t>P2,t) plane would test higher-order terms of the standard collinear pQCD. These are 
also regions where UGDFs can be tested, provided that not too big transverse momenta of 
jets taken into the correlation in order to assure the dominance of gluon-initiated processes 
(for larger transverse momenta and/or forward/backward rapidities one has to include also 
quark/ ant iquark initiated processes via unintegrated quark/ ant iquark distributions). 

Can we gain a new information correlating the space of azimuthal angle and the space 
spanned by the lengths of transverse momenta (pi t t,P2,t) ? I n particular, it is interesting how 
the jet azimuthal correlations depend on a region of (pi J t,P2,t)- F° r this purpose in FigJTBl 
we define several regions in (pi,t,P2,t), called windows, for easy reference in the following. 
They have been named for future easy notation. In FiglTJ] we show angular azimuthal 
correlations for each of these regions separately. While at small transverse momenta the cross 
section obtained with 2^2 fct-factorization and 2 — > 3 collinear-factorization approaches 
are of similar order, at larger transverse momenta and far from the diagonal pi t t = p2,t 
the cross section is dominated by the genuine next-to-leading order processes. In these 
regions the standard higher-order collinear-factorization approach seems to be the best, and 
probably the only, method to study dijet azimuthal-angle correlations. 

Cuts on pi it and p^t remove a big part of soft singularities, leaving only region of pi tt ~ p 2 ,t • 
In order to eliminate the regions where the pQCD calculation does not apply we suggest to 
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FIG. 9: Two-jet correlations for a 2 —* 3 g<? — ► component for RHIC energy W = 200 GeV. The 
soft singularities are shown as ridges. The pQCD calculations are reliable outside of the regions of 
ridges. 

exclude the region shown in diagram [15] which is equivalent to including the following cuts 
on the lengths of transverse momenta of the jets taken into account in the correlations: 

\p 1: t-p2,t\ > A s (3.1) 

In FigJT6]we show the distribution of the cross section in azimuthal angle for different (scalar) 
cuts A s = 0,2,5 GeV. We have also tried another way to remove singularities: 

\pi,t + P2,t\ > A, (3.2) 

In Fig{17|we show the distribution of the cross section in azimuthal angle for different (vector) 
cuts A„ = 0,2,5 GeV. These results are very similar to those obtained with scalar cuts. 

Both scalar and vector cuts remove efficiently the singularity of the collinear 2 — ^3 con- 
tribution at cf)_ = n. If too big values of A s or A v are used the cross section of the 
fc t -factorization 2— >2 contribution is reduced considerably. 



IV. DISCUSSION AND CONCLUSIONS 

Motivated by the recent experimental results of hadron-hadron correlations at RHIC 
we have discussed dijet correlations in proton-proton collisions. We have considered and 
compared results obtained with collinear next-to-leading order approach and leading-order 
fc r factorization approach. 

In comparison to recent works in the framework of fc t -factorization approach, we have 
included two new mechanisms based on gq — > gq and qg — > qg hard subprocesses. This was 
done based on the Kwiecihski unintegrated parton distributions. We find that the new terms 
give significant contribution at RHIC energies. In general, the results of the fc r factorization 
approach depend on UGDFs/UPDFs used, i.e. on approximation and assumptions made in 
their derivation. 
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An interesting observation has been made for azimuthal angle correlations. At relatively 
small transverse momenta (pt ~ 5-10 GeV) the 2 — > 2 subprocesses, not contributing to the 
correlation function in the collinear approach, dominate over 2 — > 3 components. The latter 
dominate only at larger transverse momenta, i.e. in the traditional jet region. 

The results obtained in the standard NLO approach depend significantly whether we 
consider correlations of any jets or correlations of only leading jets. In the NLO approach 
one obtains -W- = if 0_ < |7r for leading jets as a result of a kinematical constraint. 

Similarly dpi d ° dp2 t = if p ltt > 2p 2 , t or p 2 , t > 2p ht . 

There is no such a constraint in the fc r factorization approach which gives a non- 
vanishing cross section at small relative azimuthal angles between leading jets and 
transverse-momentum asymmetric configurations. We conclude that in these regions the 
/^-factorization approach is a good and efficient tool for the description of leading-jet corre- 
lations. Rather different results are obtained with different UGDFs which opens a possibility 
to verify them experimentally. Alternatively, the NLO-forbidden configurations can be de- 
scribed only by higher-order (NNLO and higher-order) terms. We do not need to mention 
that this is a rather difficult and technically involved computation. 

On the contrary, in the case of correlations of any unrestricted jets (all possible dijet com- 
binations) the NLO cross section exceeds the cross section obtained in the fc r factorization 
approach with different UGDFs. This is therefore a domain of the standard fixed-order 
pQCD. We recommend such an analysis as an alternative to study leading-jet correlations. 
In principle, such an analysis could be done for the already collected Tevatron data. 

What are consequences for particle-particle correlations measured recently at RHIC re- 
quires a separate dedicated analysis. Here the so-called leading particles may come both 
from leading and non-leading jets. This requires taking into account the jet fragmentation 
process. We leave this analysis for a separate study. 



V. APPENDICES 



A. Matrix elements for 2^2 processes with initial off-shell gluons 

In this paper we shall include the following 2^2 processes with at least one gluon in 
the initial state: 

(a) 99 -> 99, (b) gg -> qq, (c) gq -> gq, (d) qg -> qg, i.e. processes giving significant 
contributions for inclusive jet production at relatively small jet transverse momenta and 
midrapidities [3]. The last two processes were not included in Refs. (l7| . 24]. We shall show 
that at RHIC energies they give contributions similar (or even larger) to the contribution of 
the asymptotically dominant gg — > gg subprocess. 
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The matrix elements for on-shell initial gluons/partons read (see e.g. [20]) 

9 4 / iu su si 
2 9s V 3 ~ P ~ ~£ ~ 

1 , ( iu Au At i u 
8 9 H 6 F + 37 + 3fi +3 I + 3 i 



i yy yy 


2 _ 




2 _ 


l-M-gq^gq 


2 _ 


qg-^qg 


2 _ 



^'-9^^ + i 2 

For on-shell initial gluons (partons) s + t + u = 0. 

The matrix elements for off-shell initial gluons are obtained by using the same formulas 
but with s, t, u calculated including off-shell initial kinematics. In this case s+t+u = kf + k^, 
where k\ , k 2 < are virtualities of the initial gluons. Our prescription can be treated as 
a smooth analytic continuation of the on-shell formula off mass shell. With our choice of 
initial gluon four-momenta k\ = —k\ t and k\ = —k\ t . 

In Refs. jl7, 22| another formula which includes off-shellness of initial gluons was presented 

da f d 2 k^ t d 2 k 2)t 2 da 2 

T{x x ,k n )— J r {x 2j k 2t ) , (5.2) 



As 2 + u 2 u 2 + s 2 



d 2 pi, t d 2 p2,tdyidy 2 J tt n ' d 2 p u d 2 p 2 

where 2 

d^p£k^ = 2 W^T) alM 14M 52{kt + kt " Pht " M A • (5 - 3) 

The factor A is a function of momenta entering the hard process A = A(s, i, u, k± :t , k 2yt ) (see 
(l7|). The factor A has been rederived recently in Ref. 24] and the result of Leonidov and 
Ostrovsky was confirmed. 

Please note a different convention of UGDF in our paper (JF) with those in Refs. 17, 22| 
(/). The UGDFs in the two conventions are related to each other as 

F(x,k 2 )=f(x,k 2 )/k 2 . (5.4) 

In order to eliminate the delta function in Eq. fl5.3p we can use the same tricks as in the 
previous section. 

The formula of Leonidov and Ostrovsky is equivalent to our formula if we define: 

TV 2 — 1 A 

\M\ 2 off _ shell = 16 ^W2*) 2 ^^«. 2 ^r • (5-5) 



B. Matrix elements for 2^3 processes 



In this subsection we list the squared matrix elements averaged and summed over initial 
and final spins and colors used to calculate the contribution of the 2 — > 3 partonic processes 
(For useful reference see e.g. [20|, |2l|). 
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For the gg — > ggg process (fci + fc 2 - > k 3 + ^4 + ^5) the squared matrix element is 



\M\ 2 = -g 



2 ys N 2 C - 1 

[(12345) + (12354) + (12435) + (12453) + (12534) + (12543) + 
(13245) + (13254) + (13425) + (13524) + (12453) + (14325)] 



(5.6) 



^3) > 



i<j 



i<j 



where (ijlmn) = (kikj) (kjki)(kik m )(k m k n )(k n ki). 

It is useful to calculate matrix element for the process qq — > ggg. The squared matrix 
elements for other processes can be obtained by crossing the squared matrix element for the 
process qq -> ggg (p a + p b ^ h + k 2 + k 3 ) 



\M\ 2 = o 6 — 

3 

's aib 2 + a 2 b 1 a 2 b 3 + a 3 b 2 a 3 bi + a\b 3 



2iV 4 ^ a 3 b 3 (a 1 b 2 + a 2 &i) ^ aibi(a 2 b 3 + a 3 6 2 ) a2b 2 (a 3 b 1 + ai& 



(fefa) ' (k 3 ki)(kik 2j 

where the quantities tjj and 6j are defined as: 



(hk 2 )(k 2 k 3 ) 



h = ivbh 



The matrix element for the process gg 
appropriate crossing: 

9 

**, / '- **) ^7 • - w 



995 



\M\ 2 



—ki, —k 3 , —ki, —k 2 , k$ 



We sum over 3 final flavours (f = u, d, s) 
For the qg —>■ qgg process 



\M\ 2 qg _> qgg (k h k 2 ,k 3 ,h,k 5 ) = 
and finally for the process gq — > qgg 



\M[- 



qq^ggg 



(ki, -k 3 , -k 2 , k^ k 5 ) 



\M\ 



91^199 



{k u k 2 ,k 3 ,k^h 



qq^ggg 



-k 3 ,k 2 ,-k x ,k^k h ) . 



(5.7) 



(5.8) 

is obtained from that of qq — > ggg by 



(5.9) 



(5.10) 



(5.11) 



The squared matrix elements are used then in formula (I2.24j) . The contributions with 
two quark/antiquark initiated processes are important at extremely large rapidities. They 
will be neglected in the present analysis where we concentrate on midrapidities. 
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C. Running a, 



The treatment of the running coupling constants in 2 — > 2 and 2^3 subprocesses is 
important in numerical evaluation of the cross section. 

For the 2^2 case we shall try several prescriptions: 
(«i) a 2 = a s {p\ t )a s {p 2 t ), 



M a 2 = ct(&¥& 

>,2 _ ^,2 



(a 3 ) a 2 = a 2 s (pi,tP2,t)- 

Analogously for the 2^3 case: 
a 2 s = a s (j)\ t )a s {p\ t )a s {p\ t ), 



D. Jet separation 

In order to make reference to real situation, as in experiments, one has to take care about 
separation of jets in the azimuthal angle and rapidity space. 

In the case of fc r factorization calculation, when there are only two explicit jets we impose 
the following jet-cone condition: 

R12 = \l (A0i 2 ) 2 + (yi - V2) 2 < Ro ■ (5.12) 

Of course in this case A0 12 = 0-- -Ro is an external parameter. For reasonable values of 
R < 1 the condition may be active only for small <p_. We discuss the role of the extra cut 
in the paper. 

In the case of 2 — > 3 subprocesses one has to check two extra conditions: 

#13 = V^lS^+O/l-^) 2 < Ro , 

R23 = V / (A0 23 ) 2 + (s/2 - Vs) 2 < Ro ■ (5.13) 

Here one can expect slightly more complicated situation. Those two cuts reduce the corre- 
lation function everywhere in 0_ = A0 12 . 
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FIG. 10: Two-dimensional distributions in pu and p2t for KL (left upper), BFKL (right upper), 
Ivanov-Nikolaev (left lower) UGDFs and for the gg — > ggg (right lower). In this calculation -4 
< Vi,V2 < 4. 
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FIG. 16: Azimuthal angular correlations for different values of the parameter A s = 0, 2, 5 GeV. 
Different UGDF are used. The notation here is the same as previously. The jet-cone radius R\2 = 
1 was used in addition to separate jets. The notation here is the same as in FigJSJ 
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FIG. 17: Azimuthal angular correlations for different values of the parameter A v = 0, 2, 5 GeV. 
Different UGDF are used. The notation here is the same as previously. The jet-cone radius R\2 = 
1 was used in addition to separate jets. The notation here is the same as in FigJHJ 
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